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ABSTRACT 

F turnoff stars are important tools for studying Galactic halo substructure 
because they are plentiful, luminous, and can be easily selected by their photo- 
metric colors from large surveys such as the Sloan Digital Sky Survey (SDSS). 
We describe the absolute magnitude distribution of color-selected F turnoff stars, 
as measured from SDSS data, for eleven globular clusters in the Milky Way halo. 
We find that the Mg distribution of turnoff stars is intrinsically the same for all 
clusters studied, and is well fit by two half Gaussian functions, centered at /i = 
4.18, with a bright-side a = 0.36, and with a faint-side a = 0.76. However, the 
color errors and detection efficiencies cause the observed a of the faint-side Gaus- 
sian to change with magnitude due to contamination from redder main sequence 
stars (40% at 21st magnitude). We present a function that will correct for this 
magnitude-dependent change in selected stellar populations, when calculating 
stellar density from color-selected turnoff stars. We also present a consistent set 
of distances, ages and metallicities for eleven clusters in the SDSS Data Release 
7. We calculate a linear correction function to Padova isochrones so that they 
are consistent with SDSS globular cluster data from previous papers. We show 
that our cluster population falls along the Milky Way Age-Metalhcity Relation- 
ship (AMR), and further find that isochrones for stellar populations on the AMR 
have very similar turnoffs; increasing metallicity and decreasing age conspire to 
produce similar turnoff magnitudes and colors for all old clusters that lie on the 
AMR. 

Subject headings: Galaxy: halo — (Galaxy:) globular clusters: general — Galaxy: 
structure — methods: data analysis — stars: statistics — Surveys 
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1. Introduction 

During tlie past decade, Galactic stars from the Sloan Digital Sky Survey (SDSS) have 
been used not only to discover substructure, but also to trace the density of the smooth 
components of the Galaxy: the disks and spheroid. Because the sample of stars is large and 
the colors of stars are well measured, it has been possible to use the statistical distributions 
of star magnitudes to infer the distances to Milky Way substructure, without the need 
to know the absolute magnitude of every star. In an early application of this technique 
( iNewberg et al.l 120021 ) . an estimate of the most common absolute magnitude of very blue F 
turnoff stars in the Sagittarius dwarf spheroidal galaxy was compared with the most common 
absolute magnitude of similar stars in the tidal debris stream to determine the distance to 
the tidal stream. Even though there is a two magnitude spread in the absolute magnitudes 
of stars near the turnoff, these stars were used as distance indicators to discover stellar 
substructure and to estimate the distance to it. 



Later studies such as iBelokurov et al.l (120061 ): lJuric et al.l (|2008[ ) used photometric par- 



allax of a larger color selection of the stars in the sample, assuming they were overwhelmingly 



to discover substructures. Some searches for stellar halo substructure ( 


Rockosi et al. 


2002; 


Grillmair & Dionatos 


2006a 


b; 


Grillmair & Johnson 


2006; 


Grillmair 


2006, 


20091) have used a 



matched filter technique to discover stars that have been tidally stripped from dwarf galaxies 
and globular clusters. 



Since the iNewberg et al.l (|2002[ ) paper, Newberg & Yanny have continued to use the 



apparent magnitude of the turnoff stars in a discovered 
timate the distances to discovered overdensities of stars (lYannv et al. 



2004J: INewberg fc Yannvl l2006l: INewberg et al.l 120071 : IWillett et al.ll2009l : 



2003I: 


Freese et al. 


Yannv et al. 


2OO9I; 



Newberg et al.ll2010l ). Sometimes the assumed absolute magnitude for F turnoff stars is the 
same as was calculated for the Sagittarius dwarf galaxy, and sometimes it is adjusted slightly 
based on available information about the stellar population of the substructure being stud- 
ied. In general, the distance to the structure under study is estimated from the peak of the 
apparent magnitude distribution of the bluest main sequence stars. 



This technique was improved in ICole et al.l (|2008[ ). in which a maximum likelihood 
technique for measuring the spatial distribution of tidal debris from the Sagittarius dwarf 
galaxy was described. In this paper, the distribution of magnitudes for F turnoff stars with 
0.1 < {g — r)o < 0.3 was represented by a Gaussian with center M^^ = 4.2 and dispersion 
CTgg = 0.6. QSOs were eliminated from the sample with {u — g)o > 0.4. This distribution is 
approximately the same as the distribution of s imilar ly selected F t urnoff stars for three glob- 
ular clusters measured by INewberg fc Yanny! (120061 ) , except that INewberg fc Yanny! (12006! ) 
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fit a slightly larger tail on the faint side of the peak than the bright side of the peak. With 
the maximum likelihood technique, the star number counts expected for a particular stream 
and smooth halo model (assuming a single absolute magnitude value for all turnoff stars) 
were broadened in the Qq direction by the absolute magnitude distribution for the turnoff 
stars, before the number counts were compared with the data. Using this method, one could 
estimate the depth of the tidal debris stream along the line of sight, and could make a more 
accurate fit to the smooth distribution of stars in the stellar halo. 

Implicit in all of these studies has been the assumption that the absolute magnitude 
distribution of F turnoff stars is well represented by a single Gaussian distribution with a 
fixed peak and width. One might expect that the actual distribution might be a function of 
age and metallicity of the stellar component. In fact, the distribution might be different as 
a function of distance along a tidal debris stream, and different across the stellar halo and 
disks. To effectively use F turnoff stars as standard candles, one should study the variation 
in the distribution of F turnoff star absolute magnitudes as a function of stellar population. 

In this paper, we analyze eleven globular clusters in the Galactic halo, and find that the 
peak of the distribution of absolute magnitudes of F turnoff stars is typically Mg — 4.18, 
and that the asymmetric distribution can be approximated by a half Gaussian on the bright 
side with a width of o" = 0.36, and a half Gaussian on the faint side with a width of a = 
0.76. This distribution is surprisingly similar for all of the globular clusters studied, which 
range in age from 9.5 to 13.5 Gyr and from [Fc/H]= —1.17 to [Fc/H]= —2.30. We explain 
this surprising result by showing that it is consistent with the agc-mctallicity relationship for 
Galactic stars. Older clusters should have fainter, redder turnoffs. However, older clusters 
also contain fewer metals, which pushes the turnoff brighter and bluer. Although the color 
of the turnoff varies shghtly from cluster to cluster, the absolute magnitude of the turnoff 
only shifts about 0.1 magnitudes from the mean. 

Although we find the turnoff magnitude to be similar for the clusters studied, observa- 
tional effects may considerably change the properties of turnoff star distributions. As one 
samples stars to the limit of the survey, the photometric errors increase. Although these 
photometric errors are small compared to the uncertainty with which we know the absolute 
magnitude of a single F turnoff star, they can be large compared to the width of our color 
selection box. For bright magnitudes, the actual colors of most of the stars selected are 
actually within the color range selected. For magnitudes near the survey limit, some stars 
that should be selected are randomly measured with colors that arc too red or too blue, 
and a larger number of stars that are too red or too blue are randomly scattered into our 
selection range. The largest flux into and out of the color selection range is on the red side of 
the range, and primarily broadens the measured width on the faint side of the peak absolute 
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magnitude as one approaches the survey hmit. 

The results of this paper will make it possible to make more accurate measurements of 
the intrinsic density of F turnoff stars in the Milky Way's stellar halo. We present a function 
that can be used to estimate the observed width on the faint side of the absolute magnitude 
distribution. We also provide numbers for the fractional increase in the number counts as 
stars leak into the color box at distances where G main sequence stars are near the survey 
limit; and the fractional decrease in the number counts as F turnoff stars reach the survey 
limit and the large color errors scatter them out of the color selection range. 



2. Globular Cluster Data Selection 

Our analysis used photometric data from eleven globular clusters, taken from the SDSS 
database's seventh data release (DR7; Abazajian et al. 2009). There are a total of seventeen 
globular clusters within the SDSS DR7 footprint, but five clusters (NGC 2419, NGC 7006, 
Pal3, Pal 4, and Pal 14) were eliminated because they are too distant, and one (Pal 1) was 
eliminated because there were too few stars to obtain accurate measures of the F turnoff 
star distribution using SDSS data. A list of all clusters studied in this investigation can be 
found in Table 1. 

To determine the limits in right ascension and declination for selecting stars for each 
field, we used the SDSS SkyServer's Navigate tooljl] The image field of view was expanded 
so that the cluster was clearly visible, then expanded further until a sizable background 
distribution was also contained within this view. In general, the entire field of view spanned 
a rectangle with sides of lengths between 6 and 8 times the apparent radius of the cluster. 
The bounds of this rectangle were then used as the right ascension and declination limits in 
the Casjobs data query. 

Within these limits, we selected all objects classified as STAR, and that had {u — g)o > 
0.4. The latter re quirement was designed to avoid possible contamination from quasars 



( lYanny et al.l 120001 ). By selecting from the database of "STAR"s, we ensured that we ob- 
tained only one instance of each object. We extracted the extinction-corrected (denoted by 
the subscript '0') point-spread-function (PSF) apparent magnitudes with errors. We deter- 
mined Tcius, the radius within which the majority of the stars belong to the globular cluster, 
and Tcut, the radius outside of which there is little contamination from cluster stars, by vi- 
sual inspection of the data. Stars with rdus < r < Vcut were removed from the data set. An 



Located at 



http: // cas . sdss . org /dr 7/en/tools/ chart / navi . asp 
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example is shown in Figure 1. The rdus and rcut values used for each cluster are included in 
Table 1. 

The following clusters are listed in Table 1, but were excluded in some of our analyses: 

NGC 5053 has a v ery low Zinn fc West meta lhcity {[Fe/H]zw84: ~-2.58) which falls 
outside the range of the ICarretta fc GrattonI (119971 ) conversion scale, and is also below the 
minimur n metallicity v a lue fo r which Padova isochrones can be generated. In the newer 
work of Carretta et al. (I2OO9I), an \Fe/H] of -2.30 is found, w hich is within the range of 
the Padova isochrones (IGirardi et al.l I2OOOI : iMarigo et al.l 120081 ) . and therefore we use this 
metallicity value for our analysis. To indicate status as a potential outlier, results for NGC 
5053 are plotted using a red-dotted series in figures showing turnoff stars properties. 

M15 (NGC 7078) is a 'core-collapsed' cluster, and so has a very compact core. It 
may also have differe nt dynamics and stellar distributions than standard globular clusters 
( jHaurberg et al.ll2010l ). The SDSS photometric pipeline does not attempt to deblend crowded 
star fields, and so information on M15 in SDSS is highly biased towards stars found on the 
edges of the cluster. We do not expect this cluster to necessarily be consistent with the other 
clusters, but we include it in analysis anyway. We indicate M15's status as potential outlier 
in figures showing turnoff star properties by plotting it with a red dotted series. 

NGC 4147 contains only slightly more stars than M92, (583 versus 334) but they are 
more concentrated around the turnoff due to a lack of stars below Mg = 6.5. This lack 
of stars is due to SDSS crowded-field photometry detection efficiency problems at fainter 
magnitudes (see Section 4). We attempt to analyze NGC 4147 with other clusters later in 
our analysis, but we expect errors to be large due to the low number of data points available. 
In figures showing turnoff star properties, a blue dotted series is used for NGC 4147 to 
indicate it's low star counts and high expected errors. 

Cluster M92 (NGC 6341) fell on the edge of the SDSS DR7 footprint, and relatively few 
stars were observed in the cluster. The number of cluster stars was large enough to produce 
a fiducial fit with large bins in magnitude (0.5 magnitudes). Therefore we were able to fit 
a modified Padova isochrone to this cluster. However, as there are so few M92 stars in our 
data, especially close to the turnoff, that M92 was omitted from the turnoff analysis. 



3. Isochrone Fitting to Determine Cluster Distances 



In order to convert our observed apparent magnitudes to absolute magnitudes, we rely on 
the measured distances to each globular cluster. Since we would like to study the effects of age 



- 6 - 



and metallicity on the absolute magnitudes of the turnoff stars, we would like measurements 
of these quantities that are as accurate and as uniform as possible for our sample of globular 
clusters. In this section, we assemble the spectroscopically determined metallicities {[Fe/ H]) 
from a single group of authors; measurements were obtained from Zinn & West (1984) and 
then converted to more m odern Carretta & Gratton scale using the conversion provided in 
(ICarretta fc GrattonI 119971 ) . Using these metallicities, we then fit ages and distances to the 
clusters in a consistent fashion using Padova isochrones. 

The Padova theoretical isochrones were fit to fiducial sequences determined from data 
for eleven clusters found in SDSS. Stars in each cluster were separated into bins, then 
the mean and standard deviation {(Tg-r) of the [g — r)o distribution in each bin was found. 
Any stars in the g^ bin with a {g — r)o value beyond 2ag_r from the mean were rejected, 
and then the {g — r)o mean and crg_r of the remaining population was found. The {g — r)o 
mean and average Mg value were accepted as a point on the fiducial sequence once the 
entire bin was within 2ag-r of the current mean. Isochrones were then fit to the fiducial 
sequences using distance and age as free parameters, while metallicity was held constant at 
the spectroscopically determined value. 

In our initial attempts to use this technique to determine cluster properties, we found 
that Padova isochrones that are good fits to both the main sequence and the subgiant 
branch require unreasonably high ages (>15 Gyr) for most clusters. The lack of agreement 
between theoretical isochrones and cluster data has been explored by pre v ious a uthors. Using 
eclipsing binary stars in the Hyades open cluster, iPinsonneault et al.l ( 120031 ) showed that 
theoretical isochrones do not match true star populations; there are dis crepancies in mass, 



luminosity, temperature, and radius. In the second paper in the series, IPinsonneault et al. 



( I2OO4J ) uses Hipparcos parallax data for the Hyades cluster to further calibrate theoretical 
isochrones, finding that offsets in color indexes are sufficient to bring a theoretical isochrone 



i n line with real main sequence data. In the fourth and final paper of the series lAn et al. 



( I2OO7I ) fits isochrones to Galactic open clusters using color corrections i n (B — V ) () as a 



function of My, while using Cepheid variables as calibration points. In lAn et al.l ( 120091 ) 



updated Yale Rotating Evolutionary Code with MARCS model atmospheres were used to 
produce ugriz isochrones which were fit to main sequences of five globular clusters, producing 
ages and distances to these clusters. 



Using similar techniques, we seek to calibrate Padova ugriz is ochrones to the 



(120091 ) results. Comparing a Padova isochrone generated from the lAn et al.l (120091 ) derived 



An et al. 



age and distance for globular cluster NGC 6205 with our derived fiducial fit, we find that 
the difference between the theoretical isochrone and the data is very nearly linear along the 
main sequence (Figure 2). Therefore, we apply a linear color correction function in [g — r)o. 
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holding Mg as the independent variable, to the Padova isochrone to bring it into agreement 
with our fiducial fit. We fit only to the main sequence and subgiant branch, and ignored 
the giant branch (~ Mg > 3.5), as the linear trend was no longer valid brighter than the 
subgiant branch. We also did not fit to the lower main sequence, as our data does not extend 
to fainter absolute magnitudes. Our color- correct ion functional fit is given by the following 
equation: 



A ((7 



-0.015 * Ma + 0.089 



(1) 



This function was applied to all of the colors in the model isochrones and then these 
models were fit to our fiducial sequences derived from the data. For each cluster, we first 
chose a Padova isochrone that, modified by Equation 1, appeared to fit the fiducial sequence 
well. We then chose additional isochrones, identical except for age, which was varied about 
our "by eye" best fit in increments of 0.2 Gyrs. A Gaussian function was fit to the residuals 
of these isochrones fits, and the mean of this Gaussian was taken as the best fit age. 

Errors in our fit ages were determined through a Hessian matrix method by comparing 
the residuals in the isochrone fit. In the limit of a single fit parameter (age), the Hessian 
error method reduces to: 



a 



R\t + 2h) + R^{t- 2h) - 2R\t) 
8h^ 



(2) 



Where R'^{t) is the residual of the isochrone fit at age t, and h = 0.2 is the step size 
in the age determination method. Using Equation 2, we were able to determine the age fit 
errors for each cluster through the use of three isochrones: the isochrone of best fit age, and 
two isochrones generated at the best fit age ±0.4 Gyr. 

Padova isochrones fit using our correction function produce a consistent set of metal- 
licities, ages, and distances to our globular cluster sample, as presented in Table 1. Cluster 
color-magnitude diagrams, fiducial fits, and modified isochrone fits are shown in Figure 3. 



We compare ou r dist ance fits to dis t ances in three other sources (jPe Angeli et al.l (120051 ). 



Harris, et al.l (Il997t) and iDotter et al.l (120101)1. and compare ou r age s to other is o chron e- 
derived ages JPe Angeh et al.l J2005h . iMarm-Franch et al.l J2009h . and lPotter et all J201oh ). 
in Figure 4. Our distances appear to be in excellent agreement with other sources. Our 
ages agree to within the formal errors for each cluster, but appear to have a small linear 
systematic offset. The ages are a very close match around 13 Gyr, but are a Gyr or two 
higher for ages of 10 Gyr, so our age scale is slightly more compressed than the ages in the 
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comparison sources. 



4. Detection Efficiency for Stars in SDSS Globular Clusters 

It is clear from Figure 3 that the cluster data is incomplete at fainter absolute mag- 
nitudes, especially amongst the farther clusters: Pal 5, and NGC 4147, NGC 5024 and 
NGC 5053. This incompleteness b egins at a brighter na agnitude than is expected from the 



SDSS detection efficiency for stars (iNewberg et al.ll2002l ). The poorer detection efficiency in 



globular clusters is due to difficulty in detecting faint sources in highly crowded star fields, 
and in particular the poor perf ormance of the SDSS photometric pipeline in this regime 



(lAdelman-McCarthy et al.l 120081 ) . For sufficiently crowded fields, the code cannot deblend 



and resolve faint stars, since they are washed out by much brighter stars. 

To quantify the cluster detection efficiency, we examined nearby clusters with relatively 
complete CMDs and compared them to the farther, incomplete CMDs under the assumption 
that there are similar absolute magnitude distributions for the stars in nearby and distant 
clusters. We select stars from the CMDs in a rectangular box with bounds 3.567 < Mg^ < 
7.567, 0.0 < {g — r)o < 0.8, and a color cut of {u — g)o < 0.4, to coincide with areas of 
completeness in the near clusters (specifically 18.0 < Qq < 22.0 for NGC 6205) that are 
incomplete for the more distant clusters. We choose three nearby clusters, NGC 6205, NGC 
5904, and NGC 5272, and found that their normalized M^^, histograms in our box are, in fact, 
similar. Therefore, we average them to create a reference absolute magnitude distribution for 
cluster stars. We assume that this histogram represents the true magnitude distribution of 
globular clusters, then compare the normalized histograms of the four most distant clusters 
(NGC 4147, NGC 5024, NGC 5053, and Pal 5) to this reference. To make the normalizations 
comparable across the incompleteness of the farther clusters, we took the average difference 
of the first five bins (which showed a similar rising trend in all of our clusters) and scaled 
the entire histogram by this value, thereby scaling the clusters by matching initial trends. 

To quantify the detection efficiency in our globular clusters, we fit the ratios between 
the reference histogram and the distant cluster histograms using a parabolic function. We 
expect these globular clusters to have the same intrinsic absolute magnitude distribution as 
the bright clusters, but are missing faint stars that were lost to the crowded-field photometry. 
The histogram residuals and functional fit can be seen in Figure 5, and the resulting function 
is given empirically as: 
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1.0 if go < 20.92 

CDE{go) = { 0.0 if > 23.57 (3) 

1.0 - 0.14 ■ (^0 - 20.92)2 otherwise 

We can use this function to reconstruct incomplete cluster star distributions, and to 
model the effect of SDSS crowded-field photometry at farther distances. Errors in the pa- 
rameter fits in Equation 3 are 0.14 ± 0.001, 20.92 ± 0.007. 



5. Absolute Magnitude Distribution of F Turnoff Stars 

5.1. Fitting the Turnoff 

We are now prepared to characterize the absolute magnitude distribution of F turnoff 
stars in old stellar populations. We select stars in the F turnoff region (0.1 < (q — r)n < 0.3), 



which is the color range used in the ph otometric F turnoff de nsity searches in I Cole et al. 



( I2OO8I ). and the range originally chosen in lNewberg et al.l (120021 ) to include stars redder than 



most blue horizontal branch stars, but bluer than the turnoff of Milky Way disk stars. We 
then build a histogram in Mg from the cluster data, using a bin size of 0.2 magnitudes over 
the range 2.0 < Mg < 8.0, which minimizes potential contamination from non-cluster stars. 
We then divide each bin by the cluster detection efficiency for the qq that corresponds to 
the bin's Mg, thereby rebuilding the component of each cluster lost due to the detection 
efficiency. 

We subtract from this histogram the expected number of field stars in each magnitude 
bin, as determined from the sky area around each cluster, scaled to the sky area of cluster 
data. If the cluster lies on the edge of the SDSS footprint, or where sections of the clusters 
were unresolved by SDSS, the areas with no data were excluded from the area used in scaling 
the background bins. The r nagnitude histograni for the background was divided by the field 



stellar detection efficiency ( iNewberg et al.ll2002l ) as a function of apparent magnitude 



The top panel in Figure 6 shows the distribution of turnoff star Mg absolute magnitudes 
in globular cluster NGC 5272. We fit each SDSS cluster histogram with a 'double-sided' 
Gaussian distribution, where the standard deviation is different on each side of the mean. 
This choice provides a good fit to the data using a simple, well-known function. We have no 
theoretical motivation behind our choice of fitting function; however, this function appears 
to effectively match the form of our data without over-determining the system. The form of 
our fit function is given by: 
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Q{x; fi, ai, ar, A) = A ■ exp 



1 {x- jjf 

2 a? 



(4) 



where: 



0"i 



o-« ifMg < /i 
cr,, ifMg > fl 



When normahzed: 



^ = ^ ^ (5) 

The fit parameters A, /i, cx;, and cr^ are the amphtude, mean, left-side standard deviation, 
and right-side standard deviations, respectively. We count all bins outside of the range 
2.0 < Mg < 8.0 to be zero. 

We assumed Poisson counting errors when fitting this function to our data histograms. 
Typically, the approximation of 'square- root n' is used to simulate true Poisson errors, but 
we found that this over-emphasized low-count bins at the expense of the overall fit. We 



instead modeled our Poisson errors in the fashion of Equation 10 from iGehreld (119861 ). At 



the standard 1-cr confidence level, this equation simplifies to 6n = 1 + \/n + 1. 

A two stage process was used in fitting the parameters of the double-sided Gaussian 
function to the Mg histograms. First a Markov-Chain Monte Carlo technique was used to 
sweep the parameter space, then the best fit point was fed into a gradient descent algorithm. 
We found that this two-stage fit allowed us to avoid local minima and find the true best 
fit. A sample histogram and functional fit are provided for M3 (NGC 5272) on the top 
plot of Figure 6. Results from the fits to ten clusters can be found in Table 2. A Hessian 
matrix, calculated at the best fit parameters, was used to determine the model errors in the 
parameter fit values. This matrix was multiplied by 2 to make it equivalent to the curvature 
matrix, then inverted, so that the diagonal elements are equivalent to the squared parameter 
variances. 

To determine the turnoff color, {g — r)To, for a cluster, we fit a Gaussian profile to 
the {g — r)o histogram of all stars within 0.5 magnitudes of the /i parameter determined 
by the turnoff fits. The mean value from this Gaussian is taken as the turnoff value of the 
cluster. This definition of turnoff color is then consistent with the turnoff color of theoretical 
isochrones (the bluest point of the isochrone turnoff). Note that we do not want to choose 
the 'bluest' (minimum {g — r)^ value) point of our data, since that ignores the intrinsic color 



- 11 - 



spread of the data, and could result in a bluer measurement of the turnoff with distance as 
color errors become larger. 

5.2. F Turnoff Distribution Results 

We now examine our F turnoff distribution fit parameters as a function of the ages, 
metallicities, and distances to the clusters in our sample. Keep in mind that throughout this 
section, we are describing the turnoff magnitude in the SDSS g filter and turnoff color in 

{g-r)o. 

In Figure 7, we find no significant correlation with cluster age and our fit parameters 
II, ai, Ur-i and {g — r)TO- Intuitively, as a cluster ages // should move to higher magnitudes 
(fainter) and the turnoff color should move to smaller (redder) values of {g — r)o. The 
invariance of these parameters with age is evidence for a fairly uniform globular cluster 
sample and for the the Age-Metallicity Conspiracy discussed in Section 7. 

When these fit parameters are plotted vs. metallicity, possible relationships emerge 
(Figure 8). In general, when the metallicity of a theoretical isochrone is increased while 
holding other variables constant, the turnoff becomes more "pinched" and dimmer - that is, 
the main sequence and subgiant branch move closer together in magnitude while the turnoff 
point becomes fainter. The decrease of the ai and Gr parameters with increasing [Fe/H] 
indicates that "pinching" is occurring, however, our ii fits decrease, showing the opposite 
behavior from what is expected. This is evidence that metallicity is not the dominant factor 
in determining turnoff brightness, and will be further explored in Section 7. A reddening 
(increasing) of turnoff color is expected with increasing [Fe/H], and is observed in our fits. 
Therefore, metallicity is dominant over age in determining the turnoff color. It is important 
to note that, for our cluster sample, the absolute magnitude distribution has little dependence 
on metallicity. 

Having explored the F turnoff distribution as a function of age and metallicity, we 
now turn our attention to the distance (Figure 9). As we do not expect position in the 
Galactic halo to significantly change the structure of a cluster, we expect all of our turnoff 
fit parameters to be intrinsically independent of the observed distance for similar clusters. 
Any parameter dependence on distance is therefore actually a result of observational errors. 
We see from Figure 9 that and ai stay roughly constant with distance, and 0"^ rises first, 
then decreases. The turnoff color {{g — r)To) appears to be invariant with distance. In the 
next section, we study the observational errors contained within the SDSS database, and 
show that these errors explain the observed changes in the distribution of F turnoff stars as 
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a function of distance. 



6. Observational Effects of SDSS Errors 

6.1. Photometric Errors in SDSS Star Colors 

Photometric errors widen the observed magnitude and color distributions of turnoff 
stars in globular clusters as plotted in color-magnitude diagrams, by an increasing amount 
with fainter magnitudes. We now examine the photometric errors in the SDSS database 
and extrapolate the observational biases that result from these errors. SDSS forgoes the 
conventional "Pogs on" logarithmic defin ition of magnitude in favor of the "asinh magnitude" 



system described in lLupton et al.l ( 119991 ). These systems are virtually identical at high signal- 
to-noise, but in the low signal-to- noise regime asinh magnitudes are well-behaved and have 
non-infinite errors. The two magnitude systems do not differ noticeably in magnitude or 
magnitude error until fainter than 24th magnitude (signal-to-noise ~2.0), so our analysis 
will be valid in both systems. Figure 10 shows the relationship between apparent magnitude 
and magnitude error in the Uq, Qq and tq passbands for stars in the Palomar 5 selection field. 
We fit an exponential function with an argument that is linear in apparent magnitude to the 
error vs. magnitude data, and find: 

e{uo) = 0.0027 + e(0-80«"-i9-2) (6) 
e((?o) = 0.00031 + e(°-^93"-20-o) (7) 



-0.000026 + e(o-80^«-i9-^) 



It is the negative constant term in the exponential that determines the magnitude at 
which the errors start to become significant. It is evident from the functional fits and Figure 
10 that, of the three studied passbands, errors in uq grow most quickly, while errors in are 
the smallest throughout. Most of our good data is brighter than 23rd magnitude, where go 
and ro magnitude errors are less than 0.3. This will have little effect on the overall appearance 
of the H-R diagrams, however, magnitude errors of this degree will have a noticeable affect 
on the observed colors, where the faint magnitude errors are on the order of the observed 
{g ~ values. 
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6.2. F Turnoff Contamination 

Since F turnoff stars arc selected tfirougfi color cuts, we expect SDSS color errors to cause 
the turnoff to be contaminated by misidentified non-turnoff stars. We took; cluster NGC 6205 
as a reference data set, as the relatively low distance (7.7 kpc) implies minimal photometric 
magnitude errors near the turnoff. Cluster NGC 6205 color errors do not become noticeable 
until Mg = 7.0 {go = 21.43); at our absolute magnitude limit of Mg = 8.0 {go = 22.43), 
the magnitude errors are close to 0.1 in go and 0.15 in ro, resulting in a maximum color 
error at Mg = 8.0 equal to 0.18. From Figure 3, we see that NGC 6205 stars between 
7.0 < Mg < 8.0 are sufficiently red that even for the maximum color error, these stars are 
statistically unlikely to to be detected in the turnoff color cut. We therefore can assume that 
NGC 6205 will be illustrative as an example of a cluster with an uncontaminated turnoff. 

In order to understand how the errors affect a cluster with increasing distance, we must 
first define a process that will allow us to view a cluster as though it were observed by 
SDSS at a farther distance, taking into account the increasing magnitude errors as apparent 
magnitudes increase. We first choose a new 'effective distance' (dcs) to observe the cluster 
at, and then perform a 'distance shift' on each star for each of uo, go and ro; the magnitude 
of each star is increased by an amount equivalent to observing the cluster at dgs instead of 
the original distance {do}- The magnitude error associated with the new magnitude is then 
derived from the appropriate error equation (Equations 6, 7 and 8). The original magnitude 
error is then subtracted in quadrature from the new magnitude error to produce the relative 
increase in error. The new magnitude value is then modified by the relative increase in 
error: we select a random value from a Gaussian distribution with a mean equivalent to the 
new magnitude, and with a standard deviation equal to the relative increase in error. This 
random value produces the new shifted magnitude, equivalent to observing the star at d^s, 
including the effects of SDSS magnitude errors. 

Having defined the distance shift process, we then separate the NGC 6205 cluster data 
into three color bins and examine bin cross-contamination as the distance increases. The 
three color bins are: the primary 'yellow' turnoff bin 0.1 < {g — r)o < 0.3, the 'red' star bin 
{{9 — ^)o > 0.3) and the 'blue' star bin {{g — r)o < 0.1). Each of these bins were treated as 
separate data sets. We then performed distance shifts on each bin at 1.0 kpc steps, up to a 
maximum deflf of 80.0 kpc. At each step we reinforced the {u — g)o color cut, then counted 
the number of stars that remained in their origin bin, and the number that had leaked into 
other color bins due to color errors. We performed this process 100 times with the NGC 
6205 cluster data, and averaged the results to smooth out potential random errors. We then 
repeated this process, but included the field stellar detection efficiency in our calculations in 
order to represent the actual observed turnoff population. 
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Figure 11 illustrates the composition of the selected turnoff stars (0.1 < {g—r)o < 0.3) as 
a fraction of the original 'true' turnoff count, and as a function of distance. This calculation 
assumed 100% detection efficiency for the shifted stars. The results of applying the field 
stellar detection efficiency is shown by the dotted lines. The trends in this figure are color- 
coded by the bin of origin; the black trend is a sum of all lower curves, and represents the 
total number of stars detected as turnoff stars for a given distance. Readily apparent from 
Figure 11 is the quick influx of 'red' stars into the turnoff between distances of (10.0-20.0 
kpc), and the constant loss of true turnoff stars with distance. 

Red G stars that lie just below the turnoff arc on the densely populated main sequence, 
and have {g — r)o values just above our turnoff color cut of 0.3. Even slight color error 
perturbations will tend to shift some these stars into the turnoff color cut, resulting in 
significant 'red' star contamination at relatively low halo distances. It is clear from Figure 
11 that this is a rapid effect that occurs at fairly low distances (10.0-20.0 kpc), but as 
the errors continue to increase it becomes just as likely for a red G star to jump over the 
turnoff selection box as it is to land in it, thereby causing the red star contamination to stop 
increasing around 20.0 kpc. 

It is interesting to note that the {u — g)o color cut causes some suppression of the red 
contamination effect. Since red stars arc fainter and farther down the main sequence, they 
have higher measured magnitudes and magnitude errors than true turnoff stars. As the 
Uq passband has the highest associated magnitude errors, faint stars with large errors are 
perturbed more in {u — g)o than in {g — r)Q. Therefore, faint red stars with large color errors 
may be perturbed beyond the {u — g)o cut and subsequently removed from the data set, 
even if errors would place that star in the {g — r)o cut for the turnoff. The {u — g)o cut then 
serves to remove some of the red contamination stars from the turnoff selection. 

As the errors increase, true F turnoff stars can only leak out of the turnoff selection box. 
Since there are very few bluer {{g — r)o < 0.1) A-type stars or redder subgiants in globular 
clusters, and these in any event are bright, few A stars leak into the F turnoff star selection 
due to errors in color. Around the distance that the red star contamination stops (~25.0 
kpc), the number of stars in the turnoff selection box is comprised of roughly 60% true F 
turnoff stars and 40% redder G star contamination. As distances increase, the fraction of 
true turnoff stars in the color selection range decreases to 50%, and the total number of stars 
selected as turnoff stars decreases. This is a significant effect that has never been accounted 
for in previous research papers. 

So that future authors may compensate for these effects, we provide analytical functions 
for the F turnoff dissipation and the red star contamination. These fits arc to the 100% 
detection efficiency case (solid curves in Figure 11), representing the effect caused by color 
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errors only. We represent the F turnoff dissipation with a 4th-order polynomial function in 
des, and fit the red contamination with a similar function of 7th-order, with the coefficients 
given by ay and ar (where a = {ao,ai,a2, ...), subscripts corresponding to order of term), 
respectively: 

ay = (1.06, -0.031, 0.00020, 2.54 x 10"^ -2.67 x 10"®) (9) 

ar = (0.016, -0.020, 0.0066, -0.00043, 1.26 x 10-^ -1.92 x 10"^ 1.47 x 10'^, -4.54 x 10"^^) 

(10) 

These functions are valid in the range of 0.0 < d^s < 80.0, that is, the range of the 
distance shifts used in the above analysis. 

6.3. Effects of Magnitude Errors on the Turnoff Absolute Magnitude 

Distribution 

In order to study the effects of SDSS errors on our measured fit parameters, we performed 
distance shifts (see previous section) on the uq, go and tq magnitudes of each cluster at varying 
des steps, to a maximum of 44.0 kpc. At every d^s, we applied the {u — g)o < 0.4 color cut 
to remove stars that would have been removed as if we had selected this data from the SDSS 
database. 

At each distance shift we took into account the background subtraction and detection 
losses. Since distance shifts must be performed on a dataset of individual stars, while the 
correction functions must be applied to histogrammed data, we performed the distance shift 
before background subtraction and detection efficiency correction could occur. In order 
to keep the background subtraction consistent with a new, effectively more distant cluster 
data set, we applied an equivalent distance shift to the background prior to binning and 
subtraction. This reproduces the effect of subtracting the background prior to the shift. We 
then divided by the cluster detection efficiency, but with the parabola center shifted with 
the cluster to the new magnitude. If we did not shift the detection efficiency function before 
dividing, it would be applied to the wrong portion of the cluster histogram, since the cluster 
has been shifted to higher magnitudes. 



Before performing the functional fit to our shifted data, we apply one of three observa- 
tional biases to our Mg histogram: Th e cluster detecti o n effic iency detailed in Section 4, the 
stellar detection efficiency described in lNewberg et al.l (120021 ) . and 100% detection efficiency 
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(in which no correction is apphed). The first bias system reveals the evolution with increas- 
ing distance of globular cluster turnoff distributions as observed in SDSS data. The second 
bias system will produce the evolution of non-cluster turnoff distributions as a function of 
distance. The final system reveals the evolution of turnoff distributions if no detection bias 
is applied, that is, if all of the stars originally detected in a cluster continue to be detected 
as the distances increase. 

All ten clusters fit in Section 5 were distance shifted to increasingly greater d^s values, 
using the methods outlined above. At each d^s, a new set of fit parameters for the observed 
absolute magnitude distribution were evaluated using the methods described above. An 
example of the histogram and fit of cluster M3 (NGC 5272), shifted to a ciefr of 44.0 kpc, is 
presented in the lower plot of Figure 6. 

We present the results of the distance shifted fits for the parameters fi, ai, and 0"^ in 
Figures 12, 13, and 14 respectively NGC 5053 and M15 (NGC 7078) are plotted with a red 
dotted series to indicate their status as expected outhers, while NGC 4147 is plotted with a 
blue dotted series to indicate that it contains few stars and therefore the fits contain large 
errors. 



6.4. Observed vs. Intrinsic Absolute Magnitude Distribution of F Turnoff 

Steirs 

It is important to note that all of the clusters studied, including suspected outliers, have 
similar /i and o"; values despite differences in distance, age and metallicity. Although we see 
differences in cr^, we have shown these are due to photometric errors, and not differences in 
the absolute magnitude distribution of turnoff stars in globular clusters. This implies that 
the halo cluster population is intrinsically similar throughout. In the next section we will 
show that this similarity is related to the Age-Metallicity Relationship. 

In Figure 15 we show a series of fits to the turnoff star magnitude distribution for 
nearby cluster NGC 6205 at increasing effective distances, including the effects of the cluster 
detection efficiency. From this plot, the most obvious effect is the loss of stars with distance; 
however, one can see that these losses balance to cause ai to stay constant throughout, and 
fi shifts slightly to the left (brighter magnitudes) as the detection efficiency cuts in with 
distance. 

We find that fi is approximately constant with distance, regardless of the applied de- 
tection efficiency bias. From the plot of fi vs d^s (Figure 12), we see that ^ values decrease 
slowly with increasing distance; however, fit errors increase quickly with distance due to 
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the loss of turnoff stars. Within the fit errors, /j, is adequately described by a constant value. 
The error-weighted constant- value fit to n gives /x = 4.18(± 0.008), with a cluster dispersion 
of 0.073. We also provide a linear fit to fj, with distance, and find that a slope of zero (which 
would imply a constant value) is within the one-cr errors: 

Kdes) = -0.011(±0.02)defr + 4.39(±0.57) (11) 

We find that ai values also stay constant with distance, regardless of the applied de- 
tection efficiency, to within calculated errors (Figure 13). It is also apparent that the two 
expected outher clusters NGC 5053 and NGC 7078 have ai behaviors that differ from the 
other clusters. An error- weighted average to the ai values, excluding the two outliers, give 
ai — 0.36(± 0.006), with a cluster dispersion of 0.18. 

The values of cr^ do not stay constant with the distance shifts for any of our three 
detection efficiency cases. Because of increasing color errors with distance, all but the nearest 
turnoff star populations are contaminated by redder main-sequence stars, as described in 
Section 6.2. These redder stars enter the turnoff histograms on the fainter (> n) side, thereby 
widening the overall distribution and increasing a,, while leaving the other two parameters 
unchanged. The nearest clusters (NGC 5272, NGC 5904, NGC 6205, NGC 7089; excluding 
the core-collapsed NGC 7078) are near enough that they do not exhibit significant red main- 
sequence contamination, and are consistent with each other. We report that these nearby, 
uncontaminated clusters are representative of the "true" globular cluster distribution, with 
a (Tr fit of 0.76(± 0.04), equivalent to the fit value of 7, below. 

When we apply the SDSS detection efficiency for globular clusters, we find that all of 
our clusters show the same behavior as a function of distance (Figure 14, left). The initial, 
quick rise in cr^ with distance is due to a large influx of red main-sequence stars due to 
color errors. If the cluster is observed at even farther distances, ar is reduced due to the 
cluster detection efficiency removing increasingly larger portions of the faint edge of the 
turnoff. This consistent series in ar is evidence that the observed spread in initial cluster 
ar values is not a real feature of the clusters, but is instead an observational bias due to 
the incompleteness of SDSS crowded-field photometry at faint magnitudes. In Figure 16, we 
show a fit to the variation of cr^ with distance. 

In order to see what happens to the distribution of of turnoff stars as a function of dis- 
tance for SDSS field stars, we also apply the SDSS stellar detection efficiency when distance 
shifting the cluster. We find that the nearby clusters (NGC 5272, NGC 5904, NGC 6205, 
NGC 7089, and NGC 5466) follow a similar pattern as in the cluster detection efficiency 
system, while the initially more distant clusters (NGC 4147, NGC 5024, NGC 5053, Pal 5) 
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maintain a relatively constant value for Gr- The initially more distant clusters are already 
incomplete due to the cluster detection efficiency, which has modified their apparent cr^ fit. 
For the nearby clusters, the initial rapid rise in with distance peaks at a higher value of a^, 
and the subsequent gradual decline is less severe. This is due to the wider and fainter drop 
off for the stellar detection efficiency as compared to the cluster detection efficiency; that is, 
the stellar detection efficiency starts removing fainter stars at greater distances, and to a less 
severe degree. As the stellar field detection efficiency is the dominant observational bias in 
the SDSS, we provide a 4th order polynomial fit to versus distance, with the coefficients 
given by Hsdss: 



asdss = (-1-7, 0.46, -0.02, 0.00057, -4.7 x 10"^) (12) 

If we assume 100% detection efficiency at all magnitudes, we find that nearby clusters 
see a quick rise in cr^ as in the previous two cases, but then level out at a constant value 
(Figure 14, right; Figure 16), while farther clusters remain constant (as in the previous case). 
As discussed in Section 6.2, the ratio of "true" turnoff stars to red contaminants remains 
constant after the initial influx of red stars into the turnoff color cut. Since this ratio remains 
constant, there is no appreciable change to the observed turnoff distribution after the initial 
rise. Figure 16 shows the difference in ar evolution for different detection efficiency cases. 

The sudden inflow of red turnoff contaminants between 10.0 and 15.0 kpc is responsible 
for the rapid rise in flts for all three detection efficiency systems. The nearby clusters 
(NGC 5272, NGC 5904, NGC 6205, and NGC 7089) are close enough (< 11.5 kpc) to not 
contain significant turnoff contamination; they are therefore representative of the intrinsic 
Or value for globular clusters. When we assume 100% detection efficiency at all magnitudes, 
the deviation of cr^ fits from this intrinsic value is then purely an effect of the color errors 
due to magnitude, and are not infiuenced by pre-existing red star contamination or detection 
efficiency losses. These ar trends then serve as a basis for understanding the effects of color 
errors on the observed distribution of turnoff stars, and so we provide a sigmoid functional 
fit to the 100% detection efficiency case: 



a 

'"^ - 1 + e-(<^ee-^) + ^ (13) 

Where rfes is the effective distance of the shifted cluster. We find best fit values for the 
sigmoid functional fit to the nearby clusters are a = 0.52(zb 0.04), /3 = 12.0(± 0.31), and 
7 = 0.76(± 0.04). The three different trends for (T^ with distance, produced by the three 
separate detection efficiency systems, are compared in Figure 16. 
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We have shown that the intrinsic turnoff distribution, as quantified by the fit parameters 
fi, ai, and cr^, is similar for all observed clusters, but that the measured value of depends 
on the distance to the cluster due to photometric errors that increase if the cluster is farther 
away. Stars that affect ai are bright, and therefore have smaller color errors than stars that 
affect (Tr, and will therefore not be shifted into or out of the turnoff color cut as easily. Also, 
there are few stars adjacent to the brighter region of the turnoff, so potential contamination 
is low. Finally, detection efficiency biases must effect the fainter turnoff stars before they can 
affect the brighter stars. Therefore, bright turnoff stars are well insulated from contamination 
and detection bias. The dominant process affecting ai, then, is the loss of turnoff stars due 
to color errors. Turnoff star losses also affect fi, but at large distances fi will also be modified 
by detection biases. 



7. The Age-Metallicity Conspiracy 



The age-metallicity relationship (AMR) describes the observed relationship between a 



cluster's age and its average metallicity. Re cent observational work (jPe Angeli et al.l 12005 



Marm-Franch et al.ll2009l : iDotter et al.ll201l[ ) indicates that there are two such AMRs present 
in the Milky Way: one in which metallicity increases with age, attributed to clusters and 
dwarf galaxies that were that were gradually captured by the Milky Way over time (consid- 
ered by these authors to be "outer halo" clusters); and a set of clusters with age 13 Gyr 
that spans a very large range of metallicities, believed to be due to old clusters that formed 
rapidly alongside the Milky Way during the initial formation event (considered to be "inner 
halo" clusters). 

In Figure 17, we show that the clusters in this study follow the trend of metallic- 
ity decreasing with age, w hich we will henceforth refer to as the AMR, as presented in 
Muratov &: Gnediru (j2010[ ). Muratov & Gnedin provide a galaxy-independent model of 
galaxy formation history using semi-analytic models that take into account cosmological 
simulations. Note that the only two clusters in our sample that are not good matches to 
Muratov & Gnedin have already been identified as outliers in previous discussions. Since all 
of our clusters are located at high galactic altitudes, and are likely members of the galactic 
halo (and though to have been accreted during ga l axy a ssembly), it is not surprising that 
our clusters are similar to the iMuratov fc Gnediru ( l2010l ) relations. Figure 17 shows that 
the Muratov & Gnedin AMR model do es not reproduce the co n stant (old, rapid co-formin g 
clusters) metallicity trend described by iDe Angeli et al.l ( 120051 ) . iMarm-Franch et al.l ( 12009! ) . 
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Dotter et al. ( 2011 lR but that it is consistent with the AMR (young, late accretion) cluster 



values from this work and other papers. 

Using the Muratov & Gnedin AMR and the modified Padova isochrones, we show that 
they predict the absolute magnitudes of turnoff stars will be similar for old stellar populations 
that follow the AMR. Taking the approx i mate mean metallicity at range of ages over 8 



Gyr from Figure 8 of iMuratov &: Gnedinl (l2010l ) (plotted in Figure 17 as the solid line) 



we produce a set of metallicity versus age points that are representative of the AMR. We 
show the Padova isochrones (modified by the linear correction function in Equation 1) in 
Figure 18. The blue turnoff point can be found between 0.20 < g — r < 0.23, the turnoff 
magnitude lies between 3.77 < Mg < 4.16, and the subgiant branch is constrained between 
magnitudes of 3.4 < Mg < 3.6. Also shown in Figure 18 are two isochrones generated for 
age-metallicity combinations far from the AMR. These illustrate that the tight grouping of 
the series is a property of clusters on AMR, and is not true for arbitrary combinations of 
age and metallicity. 

The turnoff parameter fits to the unmodified cluster data in Table 2 are consistent with 
the trends seen in Figure 18. Along the AMR, as age decreases and metallicity increases, 
we see the turnoff move to slightly brighter magnitudes (lower /i) and slightly redder colors 
(higher {g — r)o). We find that the AMR produces clusters with similar isochrones, but 
the age of a cluster slightly dominates over the metallicity in determining the go turnoff 
brightness, while the metallicity slightly dominates over the age in determining the {g — r)o 
turnoff color. 

This finding implies that age and metallicity values along the AMR "conspire" to pro- 
duce a similar population distribution for all older (8.0+ Gyr) clusters. As a cluster ages, 
it's turnoff becomes redder since it has had more time for stars to evolve off of the main- 
sequence. However, an increase in metal content also corresponds to a redder value for the 
entire isochrone. The plot of isochrones along the theoretical AMR in Figure 18 implies 
that these two effects almost exactly cancel each other out, and therefore old stars formed in 
accordance with the AMR should fall in a very narrow range on a color-magnitude diagram. 
Provided one does not select a cluster from the old ("inner halo") population, this finding 
simplifies distance measurements for the predominantly old stellar population in the Galac- 
tic halo, but complicates the age determination process, which depends on the uniqueness 
of isochrones. In this paper, we have removed this complication by using only metallicities 
determined from spectra; only age and distance were used free parameters in our isochrone 



Potter et al.l |201ol) used Zinn & West metallicities; these were converted to the Carretta & Gratton 



scale in Figure 17 to be consistent. 
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fits. 



8. Application to SDSS Data 

In ttiis paper we tiave described tlie SDSS pliotometric errors and the resulting effects 
on observed F turnoff distributions (0.1 < {g — r)o < 0.3). In this section we will provide 
a method by which SDSS observations can be corrected for these effects. If we sum the 
polynomials whose coefficients are presented in Equation 9 and Equation 10, we will produce 
a function that describes the ratio of stars selected as turnoff stars to the number of actual 
turnoff stars present, as a function of distance 

eir) = ^ = j2Mr) + ^r,ir)y (14) 

where r is the distance to the stellar population. Note that this equation assumes a (m — g)o 
cut to the data; if this cut is not performed, then additional stars will be included in the 
data set and this equation will not be valid. Instead, the following equation should be used: 

^(r) = ^ = E(byi + bH)r^ (15) 

byi = (1.02, -0.011, -0.00043, 9.7 x 10"^ -5.5 x 10"*^) (16) 
bri = (0.034, -0.039, 0.011, -0.00071, 2.0 x 10"^ -2.9 x 10"^ 2.1 x 10"^ -6.4 x lO"^^) (17) 

For simple density searches that ignore F turnoff distribution statistics, dividing the 
observed density at a given distance by e{r) is sufficient to correct data for missing turnoff 
stars due to magnitude errors. For example, if one measures the number of background- 
corrected F turnoff stars in a section of the Sagittarius tidal debris stream at a distance of 
45 kpc, that number should be corrected for completeness (if necessary) and then divided 
by £:(45 kpc) in order to produce the true number of turnoff stars in that field. 

For statistical models that seek to quantify F tur noff star densities in the Galactic halo. 



a more complicated correction method is required. In lCole et al.l ( l2008l ). the authors sought 
to map the Sagittarius tidal debris stream by performing a maximum likelihood fit of a 
statistical model to F turnoff stars taken from the SDSS data. They include the effect of 
the turnoff star distribution by convolving a normalized Gaussian, representing the turnoff 
star distribution (mean=4.2, a =0.6), with a halo density model in which all of the turnoff 
stars have the same absolute magnitude {Mg =4.2). A correction for the stellar detection 
efficiency is subsequently applied. 
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We seek to modify the methods of ICole et aLl (120081 ) by incorporating our resuhs into 
their analysis. In order to include the turnoff star distribution presented in this paper, the 
turnoff distribution in the density convolution must be changed from a symmetrical Gaussian 
(A^ in that paper) to the double-sided Gaussian, Q{gQ] fj,,ai,ar, A) presented in Equations 
4 and 5, with a ar given by Equation 13. Note that we use the 100% detection efficiency 
function for ar, as the stellar detection efficiency is applied subsequent to the turnoff con- 
volution. The normalization for Q must be adjusted for the fact that turnoff stars are lost 
(and a smaller number of G stars are gained) at farther distances. To account for the turnoff 
dissipation and contamination, the distribution Q should be multiplied by 6{g). Putting 
this together, we produce a convolution recipe which provides a statistical description of the 
Galactic halo turnoff population, including the effects of the double-sided Gaussian magni- 
tude distribution and t he effects of dissipation and contamination with distance. Equation 
14 in 



Cole et al.l (120081 ) is replaced by: 



dg £ {r{g, /i)) p^{g, n)g {go - g; fi, ai, ar{r{g, /i))) 



(1^ 



where p represents a density function of magnitude and solid angle (fl), Pfj.{g, fl) is the model 
turnoff density function, in which all stars are assumed to be at absolute magnitude p (4.18), 
e is Equation 14, and Q is Equation 4 with Equation 13 as the input for cr^. e and a^. are 

QQ — l-L — 10 

both functions of r, which is equal to the distance (r = 10 § ). 

Using these examples as templates, and with the information contained in the rest of 
this paper, future authors should be able to adapt our results to correct their models for 
observational biases due to turnoff magnitude errors. Future authors must be prepared to 
compensate for three mechanisms that cause turnoff star incompleteness and misidentifi- 
cation: The combined effects of dissipation a nd contamination, giv en by Equation 14; the 
SDSS stellar field detection efficiency, given in lNewberg et al.l (|2002[ ): and the intrinsic mag- 
nitude distribution of turnoff stars within the Galactic halo, as discussed in Section 6 and 
this section. 



Discussion 



9.1. Effects of Abundance Variation 



(Pritzl et al. 


2005; 


Roederer 


2009) 



Gratton et al. 2004 and McWilliam 1997) indicate that both outer halo globular clusters 
and halo field stars have similar alpha-element abundances, consistent with an average [a/Fe] 
value of 0.3. This is strong evidence that halo clusters and field stars share a similar 
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formation history, and it is therefore reasonable to assume that they share similar population 
distributions and photometric properties. These studies also find that halo clusters and stars 
are chemically distinct from dwarf spheroidal galaxies (jVenn et al.l l2004j ) and the galactic 
disk, including stars in the thick disk. While thick disk stars can be removed from SDSS 
photometric data through magnitude and color cuts, dwarf spheroidal stars may be present 
in the halo due to tidal disruptions of infalling dwarfs. If significantly different [a/Fe] affect 
the photometric properties of a population, disrupt ed dwarf spheroid al galaxies may pose a 
problem for our technique. However, the success of ICole et al.l (120081 ) in mapping segments 
the Sagittarius tidal stream with a similar technique is evidence that that dwarf galaxy 
turnoff star distributions are not significantly different from those in the globular clusters in 
this study. 

The set of globular clusters in this study are limited to old, metal poor galactic halo 
clusters, with a maximum metallicity of [Fe/H]= —1.17 and a minimum age of 9.5 Gyr 
(NGC 5904). Figure 18 indicates that our results can be extended to an age of at least 8.0 
Gyr. From Figure 17 we see that our entire cluster sample falls within the "metal poor" 
(blue) region of the Muratov & Gnedin AMR. While our results are well-suited to our goal 
of describing the old, metal-poor halo of the Milky Way, we have not tested how far it can be 
extrapolated to stellar populations less than 8 Gyr old or more metal-rich than NGC 5904. 

The results of this paper may not apply to globular clusters with age and metallicity 
values that fall far from the Milky Way AMR. Our results may not apply to stellar pop- 
ulations in other galaxies, which have differing assembly histories and potentially different 
age-met allicity relationships. Indeed, the two unusual globular clusters that fall outside of 
the AMR in Figure 17, are the two clusters that differ from each other in turnoff magnitude 
(Figure 11). 



Recently, multiple star formation perio ds have been de tected in glob ular clusters (IMilone et al 



2008; 


Bedin et al. 


2004; 


Piotto et al. 


2007) 



multiple stellar populations to significantly affect our results, however, certain clusters with 
several separate, strongly visible main sequences (such as M54) may be poorly described by 
our results. 

Above average helium enrichment has been proposed as a poten tial mechanism for 
producing separat i ons in multiple-pop ulation cluster main sequences ( jPiotto et al.l l2007l ; 



Milone et al.ll2008l : iPasquini et al.l 120111 ). Due to the difficulty of measuring He enrichment, 
the complete effects of enrichment are still under investigation (see preliminary work by Val- 
carce at al. 2011, and Milone & Merlo 2008), but it is thought that it will result in slightly 
brighter subgiant branches and slightly bluer turnoffs, since He is less opaque than H. If the 
He abundances have only a small effects on photometry, they will not significantly influence 
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our results. 



9.2. Potential Influence of Binaries 



Binary stars present in globular clusters have the potential to bias color-mag nitude 
diagrams towards brighter, redder values, as detailed in ( iRomani fc Weinberglll99ll ). The 
maximum increase in brightness occurs when the binaries are of the same spectroscopic type, 
2.5 * log{2) ^ 0.75 magnitude, and the maximum increase in "reddening," which serves to 
widen the color-magnitude diagram towards redder values, can shift colors by as much as 
0.05 magnitudes. Depending on the binary fraction (/), this effect could introduce significant 
biases when determining population statistics of a globular cluster. 

The binary fraction in old globular clusters is generally low, arou nd a few percent to 
20%, with the bina ries being concentrated near the cluster's center ( ISoUima et al.l 120071 : 
Fregeau et al.l 120091 ) . although a select few globula r clusters have f valu es as high as 50%. 
The primary result of Monte-Carlo simulations in iFregeau et al.l ( 120091 ) is that "True" bi- 
nary fractions are likely to remain const ant with age, i mplyi ng that we can also rule out 
age-dependent biases on cluster / values. ICarney et al.l (|2003[ ) studied the binary fractions 
of metal-poor ([Fe/H] < -1.4) field red giants and dwarfs, and found them to be 16%±4% 
and 17%±2%, respectively, which is on the high end of a typical globular cluster /. Glob- 
ular clusters also feature apparent binaries, which are stars that appear as binaries due to 
crowding, which become more likely as one looks closer to the densely-packed cluster cores. 

Due to the poor performance of the SDSS photometric pipeline in crowded fields, the 
centers of globular clusters are absent from the database, and therefore ar e not present 
i n thi s study. If binaries a re concentrated towards a cluster's center as in ISoUima et al. 
(120071 ): iFregeau et al.l (120091 ). then this "disadvantage" serves to reduce the potential bias 
that binaries would have on our results, and will also greatly reduce the number of apparent 
binaries in our data. 

The efi^ect of a binary population is to shift the observed turnoff to brighter magnitudes, 
and to increase the spread in magnitudes. We ran an experiment in which we generated 
1000 stars with magnitudes given by a double-sided Gaussian with fi = 4.2, ct; = 0.36 
and cTj. = 0.76, and with colors given by a standard Gaussian with a /i(g-r) = 0.25 and 
o'(g-r) = 0.025. We then randomly selected 20% of these stars and shifted them brighter by 
0.75 magnitudes and redder by 0.05 magnitudes, which should produce a larger effect than 
we would expect from the binary fraction of a typical halo globular cluster. We found that 
the 20% binary population had a new magnitude fit of /i = 4.22, o"/ = 0.55 and cr^ = 0.75, 
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and a color fit of iJ.(g-.r) = 0.26 and cr{g-r) = 0.03. Tlie variations in tliese values are all within 
the errors of the fitting process, except for ai in the magnitude fit. It is possible that the 
quantity o"/ is sensitive to to the binary fractions typically found in halo globular clusters, in 
the sense that a large binary fraction could produce a larger cx;. The spread in o"/ values in 
our accepted sample is low (Figure 13), so we do not expect that the effect of binaries has 
changed our results significantly. A future study could potentially use a more rigorous study 
of binary fraction effects to see if the binary fraction is directly correlated to the value of ai. 

We expect the results of this paper to provide a close approximation of the F turnoff 
absolute magnitude distribution in old, metal poor stars in the Milky Way halo, as observed 
by SDSS. These results should also be a reasonable approximation of small dwarf galaxy 
turnoff distributions. We do not expect the results to be applicable to young, relatively 
metal-rich stellar populations, to populations that are outhers on the AMR, or to populations 
that are outside of the Milky Way. 



10. Conclusions 



In this paper we analyze SDSS data for 11 globular clusters in the halo of the Milky Way, 
and draw conclusions about the intrinsic and observed absolute magnitude distributions of 
F turnoff stars. The major conclusions are: 

(1) The completeness as a function of magnitude in SDSS stellar data is different in 
crowded fields such as globular clusters than in typical star fields. This incompleteness begins 
at brighter r nagnitudes and f alls off more steeply than the standard stellar detection efficiency 
reported in iNewberg et al.l ( 12002| ). We calculated an approximate parabolic function to 
describe this crowded-field incompleteness, as presented in Equation 3. 

(2) We calculated a linear correction to the turnoff region of Padova isochrones in SDSS 
colors in order to place th em on an age scal e that is consistent with determinations from other 
authors. The findings of lAn et al.l (|2009[ ) were used as a basis for this function. Without 
this correction, Padova ugriz isochrones imply ages greater than the current measured age 
of the universe. We find the appropriate correction to be A{g — r) = —0.015 * Mg + 0.089, 
valid for 3.5 < Mg < 8.0. 

(3) Using our Padova isochrone correction, a uniformly determined set of metallicity 
{[Fe/H]), distance (kpc), and age (Gyr) measurements were calculated for 11 globular clus- 
ters observed in the Sloan Digital Sky Survey. These results are collected in Table 1. 



(4) As color errors increase with distance, the fraction of F turnoff stars detected in a 
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narrow color range decreases rapidly, while redder star contamination becomes significant. 
Using Equations 9 and 10, future authors will be able compensate for both effects. 

(5) Across a range of old stellar populations (—2.30 < [Fe/H] < —1.17, dex; 13.5 < 
Age < 9.5, Gyr), the distribution of SDSS g filter absolute magnitudes of the turnoff stars 
(0.1 < {g — r)o < 0.3) is approximately constant: jjt, — Mg — 4.18± 0.008(dispersion=0.073), 
and if the distribution is fit with two half Gaussians the brighter half Gaussian has cii = 0.36± 
0.006(dispersion=0.18), and the fainter half Gaussian, when not affected by observational 
errors, has = 0.76± 0.04. 

(6) Although the Mg turnoff absolute magnitude distribution is intrinsically the same 
for old stellar populations in the Milky Way, observational errors produce biases which 
introduce a significant difference in the absolute magnitude distribution for stars selected 
in a narrow color range. While the average and bright-side half Gaussian fit parameters 
stay relatively constant to within errors, the observed fainter-side half Gaussian width (cr^) 
changes dramatically with increasing distance (and therefore increasing errors) and depends 
strongly upon observational biases. We found that Equation 13 provides a good description 
of the change in cr^ with distance, as influenced by color errors only. This must be accounted 
for if F turnoff color cuts are to be used to trace stellar structure in the Milky Way. 

(7) The turnoff properties for old stellar populations is consistent with the known Age- 
Metallicity Relationship for the Milky Way. Older globular clusters should have fainter, 
redder turnoffs, but because they are more metal-poor (which tends to produce brighter, 
bluer turnoffs) the turnoff is similar for all older clusters. For stellar populations at least 8.0 
Gyr old, the turnoff color ranges from 0.20 < {g — r)o < 0.23, the turnoff magnitudes range 
between 3.77 < Mg < 4.16, and the subgiant branch is constrained between the magnitudes 
of 3.4 <Mg< 3.6. 

(8) We provide a method for correcting SDSS data for the effects of magnitude errors. 
We discuss two cases: (i) If the stellar population under study is at a single distance, then 
the observed F turnoff star number counts should be corrected for completeness, and then 
divided by e{r), as given by Equation 14 or 15; (ii) If the stellar population is distributed 
in distance, one can convolve the proposed density model with a function that accounts for 
observational errors before comparing with data. 

F turnoff stars have proven useful as tracers of Galactic structures, particularly in the 
halo. Before completing this work, we were concerned that the assumption of a single 
absolute magnitude distribution for turnoff stars in the halo would be problematic. Instead, 
our analysis shows that while a single turnoff is a good assumption, photometric errors 
significantly affect the stellar population as a function of apparent magnitude. Since we 
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have found that all old stellar populations are intrinsically similar, it is possible to model the 
selected populations of stars (as selected in a narrow color range) as a function of magnitude 
and correct measurements of Galactic structure for this effect. 
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Table 1. Results of modified Padova isochrone fits 



NGC 


Messier 


I 


b 


Metallicity 


Fit Distance 


Fit Age ± Error 
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'^cut 


# Stars in 
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O 


O 


CG97 [Fe/H] 


(kpc) 


(Gyr) 


o 


o 
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NGC 4147 




252.85 
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-1.58 


19.3 


11.2 ± 0.8 


0.07 


0.08 
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NGC 5024 
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332.96 


79.76 


-1.89 


18.7 


12.5 ± 1.3 


0.25 


0.30 
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NGC 5053 




335.70 


78.95 


-2.30^ 


18.5 


11.5 ± 1.1 


0.14 


0.17 
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NGC 5272 


M3 


42.22 


78.71 


-1.43 


10.4 


12.2 ± 1.1 


0.30 
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NGC 5466 
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73.59 


-2.14 


15.6 


13.4 ± 0.9 


0.18 


0.25 
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NGC 5904 


M5 


3.86 


46.80 


-1.17 


8.0 


9.4 ± 0.1 


0.21 


0.28 


3488 


NGC 6205 


M13 


59.01 


40.91 


-1.42 


7.7 


12.7 ± 0.2 


0.24 


0.33 


2192 


NGC 6341 


M92 


68.34 


34.86 


-2.17 


8.7 


13.4 ± 0.9 


0.135 


0.20 


334 


NGC 7078 


M15 


65.01 


-27.31 


-2.04 


11.0 


11.2 ± 0.9 


0.19 


0.25 


2939 


NGC 7089 


M2 


53.37 


-35.77 


-1.38 


11.5 


12.6 ± 0.3 


0.14 


0.18 


1451 


Pal 5 




0.85 


45.86 


-1.24 


21.0 


11.5 ± 0.7 


0.087 


0.12 


1478 



^The number of stars represents the number of stars contained within r^j^j, and therefore the number of stars used to build 
the cluster dataset. 

■^The Zinn &; West metallicity of NGC 5053 (-2.58) was outside the effective co nversion range of the CG97 metallicity scale, 
as well as the range of available Padova isochrones, so the lCarretta et ah value of -2.30 was used instead 



Table 2. Fit values and errors to cluster data 
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Fig. 1. — Positions of stars near globular cluster NGC 5053 in right ascension and declination. 
Stars assigned to the cluster are marked with circles, while stars assigned to the background 
are marked with crosses. At the edge of the cluster, it is difficult to separate cluster stars 
from background stars, so a ring of stars at the interface has been removed, leaving a gap on 
the plot. The empty area in the center of the plot is where the SDSS photometric pipeline 
was unable to separate individual stars, and the data was not included in the database. This 
area has been subtracted from the cluster area determinations in Section 2. 
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Fig. 2. — Diagram of Padova Isochrone fit for cl uster NGC 6205 . The values of the raw 
Padova Isochrone, (blue circles) generated from the lAn et al.l (120091 ) results, were subtracted 
from the values of the cluster fiducial fit (black diamonds) to produce the offset (crosses). A 
red vertical line indicates zero offset, or a perfect match between isochrone and fiducial fit. A 
linear fit (dotted line) to the offset was made for values below the giant branch {Mg > 3.5): 
A{g — r)o = —0.015 * Mg + 0.089. This correction causes Padova isochrones to fit the turnoff 
on color-magnitude diagrams. We then used this fit to match isochrones to our cluster 
sample (see text). 



- 3 - 




Fig. 3. — Color-magnitude diagrams, fiducial fits, and isochrone fits for all eleven globular 
clusters used in this study. Black dots represent the positions of individual stars in M^, 
{g — r)o; red circles represent 2-a rejection fiducial fits to the stars in Mg strips; and blue 
lines represent Padova isochrones modified by our correction function. The ages, metallicities, 
and distances for the isochrones are given in Table 1. Note the absence of faint stars in some 
clusters. This is due to the poor performance of the SDSS photometric pipeline in crowded 
fields. The shape of this incompleteness is described in the text. Also note the sparse data 
in NGC 6341; it will not be used for F turnoff analysis. 
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Fig. 4. — Comparison of the distance and age values from this paper with values from other 
sources. The top panel plots the distances from this paper, subtracted from the respective 
external values, plotted vs. our fit ages. Dotted trend lines are linear fits to the points. Our 
distances are in good agreement with the previous literature. The lower panel is the same as 
the top panel, but with the age differences as the y-axis. Our ages agree with the previous 
literature to within the formal errors, but have a small linear offset. Note that Palomar 5 is 
absent from this plot, as the other authors did not study this cluster. 
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Fig. 5. — The SDSS crowded-field photometry detection efficiency. Due to the poor per- 
formance of the SDSS photometric pipehne in crowded fields, more distant clusters have a 
deficit of stars at fainter magnitudes relative to brighter clusters. The ratio between the 
number of stars observed in distant clusters and the reference histogram are plotted as col- 
ored lines. The reference histogram was created from the average histogram of Mg for the 
nearby clusters NGC 6205, NGC 5904, and NGC 5272, and then shifted to the distance of 
each of the distant clusters NGC 4147, NGC 5024, NGC 5053, and Pal 5. The parabohc 
functional fit {CDE{go) = 1.0 - 0.14 ■ (^o - 20.92)^, fo r 20.92 < gn < 23.57) is plotted as the 
solid black line. The stellar detection efficiency from iNewberg et al.l ( 120021 ) is plotted as a 
dotted series for reference. 
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Fig. 6. — Histogram fits for NGC 5272 original data (top) and convolved to to an effective 
distance of 44.0 kpc (bottom) with the cluster detection efficiency applied. Blue lines outline 
the bins, while blue points mark the bin centers, to which a double-sided Gaussian (green 
line) was fit. Fit values are {ii,ai,ar) = (4.33, 0.39, 0.77) for the NGC 5272 raw data, and 
(3.97, 0.35, 0.55) if the cluster was instead observed at a distance of 44.0 kpc. Note the 
difference in total counts between the two plots; this is caused by losses due to color errors 
with effective distance. 
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Fig. 7. — Plot of globular cluster ai, ar and turnoff color values vs. isochrone fit ages. 
There are no statistically significant trends in our parameters as a function of cluster age. 
It is expected that the turnoff should shift to redder color and fainter magnitude as it ages, 
but since these clusters fall along the Age-Metallicity Relationship, the change in metallicity 
counterbalances the effects of the change in age. 
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Fig. 8. — Plot of globular cluster /i, ai, Or and turnoff color values vs. CG97 scale metal- 
licities. The fit parameters and cx; appear to decrease with metallicity, indicating that 
isochrone fits should have slightly brighter and sharply curved turnoffs as [Fe/H] increases. 
The fit parameter cr^ does not appear to be correlated with metallicity; We will show that 
changes in Ur are primarily due to distance (therefore, error) effects. One would expect an 
increase in metallicity to produce a fainter turnoff, however, since the more metal-rich clus- 
ters are also younger, the turnoff is actually slightly brighter for a higher metallicity star. A 
slight increase in the turnoff color with metallicity is observed. 
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Fig. 9. — Plot of globular cluster /j,, ai, ar and turnoff color values vs. cluster distances. 
We see that /i appears to increase (become fainter) slightly with distance, ai stays roughly 
constant, and cr^ first increases, then decreases. The turnoff color {{g — r)To) appears to 
be invariant with distance. Since cluster properties should not be dependent on a cluster's 
distance from the Sun, any trends in this figure are a result of the increasing magnitude 
errors with distance. 
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Fig. 10. — Plot of SDSS color magnitudes versus respective color errors, go and and Mq errors 
are offset by constants of 0.5 and 1.0, respectively, for illustration. Errors and dereddened 
photometric values are from the Sloan Digital Sky Survey DR7 database using stars near 
the globular cluster Palomar 5. It is clear from the plots that the errors in uq rise the 
fastest, followed by the errors in tq. uq data was cut for uq > 23.5, where the Uq errors 
became unreliable. Black lines are functional fits to the data (colored points): e(uo) = 
2.71e - 03 + eO-80«o-i9-2; e{go) = S.lle - 04 + eO-™^'o-20.o. ^^^^^ _ _2.63e - 05 + e0-80'^«-i9-8. 



- 11 - 



So 

19.2 20.7 21.6 22.2 22.7 23.1 23.4 23.7 
1.2 1 1 1 1 1 1 1 1 1 




d^ff, kpc 

Fig. 11. — Percentage of stars by type detected as F turnoff stars due to SDSS color errors 
as a function of distance/mean apparent magnitude. The yellow (middle) curve indicates 
the percentage of 'true,' (original) F turnoff stars detected as F turnoff stars, while the red 
(lower) curve indicates the percentage of 'red' stars {{g — r)o > 0.3) incorrectly identified 
as F turnoff stars due to color errors. The blue curve, which is nearly coincident with the 
X-axis, represents the same trend but for 'blue' stars {{g — r)o < 0.1). The black (top) curve 
is the sum of all lower curves, and illustrates the total percentage of stars detected as turnoff 
stars relative to the true turnoff distribution. The dotted curves indicate the effect of the 
field star detection efficiency on star counts. It is clear that the relatively higher color errors 
at fainter magnitudes, which largely consist of 'red' stars, initially causes a large number 
of red stars to contaminate the turnoff, which then tapers off as the errors increase. True 
turnoff stars quickly and continually leak out of the turnoff bin. By as near as 25.0 kpc from 
the Sun, half of all detected turnoff stars are actually redder star contaminants. 
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Fig. 12. — fi series fits to ten globular clusters, convolved to errors consistent with observing 
them at larger distances [des). The cluster detection efficiency (Equation 2) was applied 
during the distance shifts. Although the fit fi values decrease with distance, the errors 
increase as well. To within errors, the n values are consistent with a constant value. The red 
dotted series have been rejected as outliers. The blue dotted series for NGC 4147 indicates 
large expected errors due to low star counts. The error-weighted average of initial /i points 
for the remaining globular clusters is 4.18(± 0.008), plotted as the black dot-dash line. 
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Fig. 13. — CT/ series fits to ten globular clusters, convolved to errors consistent with observing 
them at larger distances {des)- The cluster detection efficiency (Equation 2) was applied 
during the distance shifts. To within errors, the o"/ values stay constant. The red dotted 
series have been rejected as outliers. The blue dotted series for NGC 4147 indicates that large 
errors are expected due to low star counts. The error-weighted average (ignoring outliers) 
of initial ai points for the remaining clusters is 0.36(± 0.006), plotted as the black dot-dash 
line. 



- 14 - 



17.8 19.3 20.2 20.8 21.3 21.7 22.0 22.3 22.6 17.8 19.3 20.2 20.8 21.3 21.7 22.0 22.3 22.6 




Fig. 14. — ar series fits to ten globular clusters, convolved to errors consistent with observing 
them at larger distances (des). The left plot represents clusters convolved while applying 
the SDSS crowded-field photometry detection efficiency function (Section 4). When this 
detection efficiency is applied, all clusters in our analysis follow the same trend (reproducing 
the observed dependence of with distance) as they are convolved to farther d^s- The right 
plot contains only clusters not initially affected by the detection efficiency, and did not have it 
applied during the convolutions, since it was not possible to remove the crowding effect from 
the more distance cluster data. This represents the ideal observation, in which no SDSS 
detection losses occur beyond the initial status of the cluster, but with only photometric 
errors increasing with distance. In both cases a quick rise in ar with distance is observed, 
due to the large influx of redder main-sequence stars into the F turnoff cut, as described 
in Figure 5. In the left plot, the detection efficiency continually cuts into the fainter stars, 
causing a constant fall in cr^ after the quick rise. When there is no loss due to detection 
efficiency (right), ar remains constant after the initial rise. To this trend we fit a sigmoid 
function (dashed line, see text): a, (3, and 7 values of 0.52(± 0.04), 12.0(± 0.31), and 0.76(± 
0.04), respectively. Dotted lines indicate clusters that have been treated as outliers. 
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Fig. 15. — Double-sided Gaussian fits to NGC 6205 turnoff (0.1 < {g — r)^ < 0.3) histograms 
with increasing effective distance (c?eff)- The third tallest Gaussian is at the original distance 
of 7.7 kpc, while the remaining Gaussians are fits to the data shifted to d^ss of 8.0 kpc 
through 28.0 kpc, in 2.0 kpc increments. Initially, the Gaussian fits become taller and 
becomes larger as redder contamination stars enter the data set. Then the peak of the 
histogram becomes smaller as the cluster moves to larger d^^ values, where increasing color 
errors cause a bleed-off of turnoff stars. Note that as the cluster is shifted beyond 20.0 
kpc (fourth curve from the bottom), the Ur value eventually begins to shrink as the cluster 
detection efficiency removes fainter stars from the data. 
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Fig. 16. — Weighted averages of series fits witli effective distance (liefr) for three different 
detection efficiency functions: no detection ef ficiency apphed , (100% detection, upper curve) 
the sigmoidal SDSS detection efficiency from iNewberg et aLl (120021 ) . (middle curve) and the 
SDSS crowded-field photometry parabohc cluster detection efficiency discussed in Section 4 
(lower curve). 
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Fig. 17. — Plot of globular clust e r age s versus metallicity, from four different sourc e s: thi s 
paper (blue circles) . IPotter et aD koid ) ; be Angeli et al.l boO^ ; iMarm-Franch et al.l feoooj ) . 
Two clusters, NGC 5053 and NGC 7078, are shown as blue rings to indicate that they are 
outliers in our analysis. T hese are overlaid on a theoretical age-metallicity relationship from 
Muratov fc Gnediru ( 120101 ) (blue- and red-shaded areas). Note that most of the clusters, 
including all of the clusters in this study besides noted outliers NGC 5053 and NGC 7078, 
are consistent with the theoretical age- metallicity relationship. N ote that there ar e some 
high age, high metallicity clusters in the iMarm-Franch et al.l (120091 ) and lDotter et al.l (120101 ) 
data that they attribute to a constant age with metallicity relationship at around 13 Gyr in 
the inner halo. 
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Fig. 18. — Plot of Padova isochrones, modified by a c orrectional fit (Equ a tion 1), from 
the Age-Metallicity Relationship (AMR) presented in iMuratov &: GnedinI (120101 ). Two 
isochrones not on the AMR are also plotted as black dotted series, at 8 Gyr and 13 Gyr, 
in order to show the behavior of extreme outliers. All of the isochrones on the AMR 
have turnoff values that are similar to each other, with a turnoff color between between 
0.2 < {g — r)o < 0.23 and a turnoff magnitude of 3.77 < Mg < 4.16. The subgiant branch 
is also well constrained in magnitude: 3.4 < Mg < 3.6. These constraints will be useful in 
determining distances to old stellar populations. The following Age (Gyr) and [Fe/H] (dex) 
value sets were used to produce the AMR isochrones: 13.0, -1.7; 12.0, -1.5; 11.0, -1.35; 10.0, 
-1.20; 9.0, -1.05; 8.0, -0.90. The age, [Fe/H] values for outlying isochrones are: 13.0, 0.0; 8.0, 
-2.0 



